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I.  INTRODUCTION 


The  transient  heat  conduction  in  semi-infinite  solids  is  an  important 
heat  transfer  problem.  Typical  examples  are  the  heating  by  propellant  gas  of 
large  caliber  gun  barrels,  impingement  heating  on  a  ship  deck  during  missile 
launching,  and  solar  heating  of  the  earth  surface.  Solutions  to  these  problems 
are  well  known  when  constant  thermal  properties^  are  assumed.  Some  materials, 
such  as  type  420  stainless  steel  and  type  4130  steel,  however,  possess  thermal 
properties  which  are  strong  functions  of  temperature  as  shown  in  Figure  1.^ 
Neglect  of  this  fact  may  lead  to  significant  errors  in  heat  transfer 
calculations.  This  is  especially  the  case  when  high  heating  rates  are 
involved. 

When  the  temperature  dependence  of  thermal  properties  is  accounted  for, 
the  heat  equation  becomes  nonlinear  and  its  exact  solution  is  unattainable. 
In  the  relatively  simple  case  in  which  the  thermal  conductivity  is 
temnerature-dependent  and  the  surface  temperature  is  constant,  Schubert,  et 
al.*^  obtained  a  similarity  solution  for  their  one-dimensional  solid  state 
model  of  the  oceanic  lithosphere  and  asthenosphere.  However,  solutions  for  a 
solid  with  both  thermal  conductivity  and  specific  heat  as  functions  of 
temperature  and  subjected  to  a  time-varying  boundary  condition  have  not  been 
reported . 

In  this  study  we  consider  a  semi- inf inite  homogeneous  solid  with  thermal 
properties  which  have  a  power  law  dependence  on  temperature.  Two  types  of 
time-dependent  boundary  conditions  are  prescribed  at  the  surface:  a  heat  flux 
“  at®  and  a  temperature  T^  ■  Tq(1  ♦  t“  ).  The  power  law  representations 
are  often  useful  in  many  engineering  applications.  Utilizing  the  method  of 
similarity  transformation  vis  one-parameter  groups,  we  transformed  the  non¬ 
linear  governing  partial  differential  equation  into  an  ordinary  differential 
equation  that  can  be  integrated  numerically  by  using  any  one  of  the  several 
ordinary  differential  equation  packages.  This  method  is  simple  and,  above  all, 
the  solution  obtained  in  terms  of  similarity  variables  provides  a  broad 
representation.  Such  an  advantage  is  usually  not  obtained  when  the  governing 
partial  differential  equation  is  directly  solved  by  numerical  methods. 

Using  type  420  stainless  steel  properties  for  sample  calculations, 
solutions  obtained  are  the  local  temperature  and  local  heat  flux  at  any 
distance  from  the  surface  of  the  medium.  In  the  limiting  case  of  constant 
temperature  or  constant  heat  flux  applied  to  the  solid  with  constant  thermal 
properties,  a  comparison  is  made  with  the  available  exact  solution.^ 
Meanwhile,  the  significance  of  the  temperature  dependence  of  thermal 
properties  for  heat  transfer  calculations  is  demonstrated. 


^H.  Carslaw  and  J.  Jaeger,  Conduction  of  Heat  in  Solids.  Oxford  Univerity 
Press,  London,  1959,  Chapter  2. 

2 

Aerospace  Structural  Metals  Handbook.  U.S.  Army  Materials  and  Mechanics 
Research  Center,  Watertown,  MA,  1980. 

'^G.  Schubert,  C.  Froidevaux,  and  D.A.  Yuen,  "Oceanic  Lithosphere  and 
Asthenosphere:  Thermal  and  Mechanical  Structure,"  J.  Geophvs.  Res.  81.  No.  20, 
1976,  pp.  3525-3540. 
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II.  FORMULATION  OF  PROBLEM 


The  equation  describing  the  transient  heat  conduction  in  a  semi-infinite 
solid  with  temperature  dependent  properties  is 

cT  g  aT 

pC  (T)  —  -  —  [K(T)  — (1) 

5t  dx  dx 


Where  T,  t,  and  z  are  the  temperature,  time,  and  distance  from  the  surface  of 
the  solid,  respectively.  Figure  2  depicts  the  solid  under  consideration.  The 
density  p  is  considered  to  he  constant.  The  thermal  conductivity  R(T}  and 
the  specific  heat  Cp(T)  are  assumed  to  be  functions  of  temperature  in  the 
following  form: 


K(T) 


o  m 


Cp(T) 


T  -  T, 


C  C  [- 
po  n*^ 
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(2) 

(3) 


where  C^,  m,  and  n  are  dimensionless  constants  which  can  be  determined 
from  the  experimental  data  of  a  given  material.  and  C  are  respectively 
the  thermal  conductivity  and  specific  heat  evaluated  at  temperature  T^.  It  is 
further  assumed  that  the  solid  is  initially  at  a  uniform  temperature  T^. 


Figure  2.  Semi-infinite  Solid  with  Temperature- 
Dependent  Thermal  Properties 


Eqs.  (2)  and  (3)  appear  unrealistic  when  T  -  since  at  that  temperature 
both  K  and  Cp  become  zero.  Practically,  however,  when  K  «  Cp  “  0,  which 
represents  the  case  in  which  the  material  does  not  conduct  heat  ^away  and  has 
no  heat  capacity,  the  temperature  will  then  jump  up  immediately  when  the 
material  is  being  heated.  Thus,  the  singular  behavior  at  T  “  appears  only 
in  a  very  short  period  of  time  and  will  have  no  significant  effect  on  the 
overall  heat  transfer  calculation  in  most  engineering  problems. 

For  the  convenience  of  analysis,  we  nondimersionalize  Eqs.  (1)  through 
(3)  by  introducing  the  following  dimensionless  variables: 
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where  L  is  the  reference  length.  Substitution  of  these  variables  into  Eq.  (1) 
leads  to  a  dimensionless  form  of  the  heat  equation 


dh  ,  as  C_  30 

0“  -  +  m0““k— e“  —  .  (5) 

?s2  as  at 


We  now  consider  two  types  of  time  dependent  boundary  conditions 
prescribed  at  the  surface  (i.e.,  at  z  «  s  «  0): 

heat  flux  *  at^  (6) 

and 

temperature  »  T^^Cl  +  t^*)  (7) 

where  a  is  a  constant  with  units  the  same  as  that  of  heat  flux  and  t  is  the 
dimensionless  time  defined  in  Equation  (4).  The  powers  b  and  b*  are 
dimensionless  constants. 

The  complete  sets  of  boundary  and  initial  conditions  for  Eq.  (5) 
corresponding  to  Equations  (6)  and  (7)  are  given  as  follows: 

A.  Prescribed  Heat  Flux  ■  at^ 

By  replacing  the  left  hand  side  of  Eq.  (6)  by 

5T 

-  -K(T)  —  (8) 

3x 


and  choosing  the  reference  length 


we  obtain  the  following  conditions  for  Eq.  (5): 


-  =  t'"  at  B  -  0,  t  >  0 
?8 

0**O  at6»®“,t>0 


(10a) 


(10b) 


0=0  atfl>0,t“0. 


(10c) 


B.  Prescribed  Surface  Temperature  ^  ^ 

In  this  case,  the  boundary  and  the  initial  conditions  are 


6=t  atc=0,t>0 

0  =  0  ats*'»,t>0 

0=0  atB>0,  t=0 


(11a) 

(lib) 

(11c) 


in  which  the  reference  length  L  can  be  chosen  to  be  any  constant  value. 

III.  METHOD  OF  SOLUTION 
A.  Prescribed  Heat  Flux  =  at^ 

We  obtained  a  solution  of  this  problem  by  using  the  similarity 
transformation  via  one-parameter  grouper”^  After  the  transformation,  the 
number  of  the  Independent  variables  was  reduced  from  two  (i.e.,  t  and  s)  to 
one  (i.e.,  V  ).  The  resulting  ordinary  differential  equation  was  then 
integrated  numerically. 


'*A.J.A.  Morgan,  "The  Reduction  by  One  of  the  Number  of  Independent  Variables 
in  Some  System  of  Partial  Differential  Equations."  Quart.  J.  Math.  (Oxford). 
2,  1952,  p.  250. 

^A.G.  Hansen,  Similarity  Analyses  of  Boundary  Value  Problems  in  Engineering. 
Prentice-Hall,  Englewood  Cliffs,  NJ,  1964. 

^W.F.  Ames,  Nonlinear  Partial  Differential  Equations  in  Engineering.  Academic 
Press,  1965,  pp.  135-141. 


In  the  transformation  procedure,  we  first  introduce  three  variables: 


t'  - 


and  6'  •  Od 


(12) 


where  d  is  a  parameter,"^,  and  are  constants  to  be  determined  by  using 
the  governing  heat  equation  together  with  the  boundary  conditions. 
Substituting  these  variables  into  Eqs.  (S)  and  (10a)  results  in 


2«2  -  (hh-Dv,  m  -  (m+Dr,  (m-1)  36'  2 
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(13) 


(14) 


To  satisfy  the  invariance  requirement  of  the  transformation,  the  indices  of 
the  parameter  d  in  each  term  of  the  equation  are  set  equal.  Then  from  Eq. 
(13)  we  have 


2«2  -  (m  +  l)/j  ■  2<»2  -  (m  ♦  D'/j  *  nj  -  (n  +  1)  Vj 


(15) 


and  from  Eq.  (14) 


"2  ”  (m  +  “  -b(»2  . 


(16) 


Solving  Eqs.  (15)  and  (16)  simultaneously  gives 

1  +  2b 

2  +  m  +  n 

“2  1  +  2b 

—  -  (m  +  1) - b. 


(17) 


2  +  m  +  n 


(18) 


and 


B  *  — 


Then  from  Eqs.  (17)  and  (18),  we  obtain 


A  -  B(m  ♦  1)  -  b. 


Now  the  similarity  variables  can  be  formulated  as 


?  -  t  ,  V 


We  seek  temperature  solutions  G  of  the  form 

yi/«i 

0(8, t)  -  t  fiv)  • 


The  temperature  function  f  in  the  equation  is  a  function  of  the  similarity 
variable  ^  only.  In  terms  of  ^  and  v,  we  obtain  the  following  expressions 
through  the  chain  rule. 


—  »  B  |®“lf  -  A  f' 

8t 

?G 

—  -  |®-Af' 

?s 

06^ 


where  the  prime  denotes  the  derivative  with  respect  to  rj.  Substituting  Eq. 
(23)  for  G  and  its  derivatives  in  Eq.  (24)  into  Eq.  (5)  leads  to  the 
following  ordinary  differential  equation. 


f"  *  f'lmf^f'  +  (Cn/C„)A’?f°"®] 

-  Bf“-®^1  -  0. 


All  of  the  constants  C  ,  C^,  m,  n,  A,  and  B  in  the  equation  are  known.  In  the 
same  manner,  the  bounoary  conditions  (10a)  and  (10b)  can  be  transformed  to 
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The  initial  condition  (10c)  is  not  needed  for  the  solution  of  Eq.  (25).  In 
fact,  the  expression  (23)  has  satisfied  the  initial  condition  (10c). 

In  the  simple  case  that  both  the  thermal  properties  of  the  solid  and  the 


boundary  heat  flux  are  constant  (i.e.,  m  *  n  *  b  ■  0, 
and  (26)  reduce  to 


1),  Eqs.  (25) 
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f  -  0 
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(28b) 

B.  Prescribed  Surface  Temperature  T^ 


Tod  +  t 


b* 


The  heat  equation  is  identical  to  Eq.  (25).  In  this  case  ve  will  choose 
a  different  set  of  expressions  for  A  and  B  defined  in  Eqs.  (19)  and  (20)  so 
that  the  boundary  conditions  (11a)  and  (11b)  are  satisfied.  Following  the 
above  procedure,  we  obtain 


and 


A  =  b*(m  -  n)/2  +  1/2 


B  =  b*. 


(29) 


(30) 
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The  boundary  conditions  for  Eq,  (25)  are 


f  -  1 
f  -  0 


at  *?  *  0 
at  ’?  ■  oo. 


(31a) 

(31b) 


The  initial  condition  is  not  needed  for  the  solution  since  the  expression 
(23)  has  satisfied  the  initial  condition  (11c). 

In  the  case  that  the  surface  temperature  T^  *  constant,  the 
dimensionless  local  temperature  6  given  in  Eq.  (4)  has  to  be  redefined  as 

T  -  T„ 
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in  order  to  yield  the  same  form  of  governing  equation  and  boundary  conditions 
given  in  Eq.  (23)  and  Eq.  (31). 


Eq.  (23)  wi'  conditions  (26)  or  (31)  is  ready  to  solve  for  the  function 
f  and  its  derivative  f'.  Once  the  values  of  f  and  f'  are  obtained,  the  local 
temperature  in  the  solid  can  he  calculated  from  Eq.  (23).  For  the  calculation 
of  the  local  heat  flux,  ve  derive  that 
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In  dimensionless  form  it  is 


(33) 
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IV.  SAMPLE  CALCULATIONS  AND  DISCUSSION 

We  chose  type  420  stainless  steel  properties  for  sample  calculations 
since  their  strong  temperature  dependence  can  serve  to  explain  the 
significance  of  the  temperature  dependence  for  heat  calculations.  Figure  1, 
reproduced  from  Ref.  2,  shows  the  experimental  data  of  the  thermal 
conductivity  R  and  the  specific  heat  C  of  the  steel.  Both  increase 
monotonically  with  temperature  in  the  temperature  range  indicated.  The  data 
can  be  nondimensionalized  through  the  use  of  the  dimensionless  variables 
defined  in  Eq.  (4).  Taking  T^  ■  294°K  (70®F)  as  the  reference  temperature, 
the  results  after  nondimensionalization  are  shown  in  Figure  3.  The 
experimental  data  can  he  approximated  by  the  dashed  lines  which  are 
represented  by 

k  -=  Cj^e®  -  00.04727  C„  -  1  and  m  -  0.04727  (35) 

Cp  -  0^6°  -  00.2491  -  1  and  n  -  0.2491 ,  (36) 


These  approximated  curves  are  much  better  representations  of  the  thermal 
properties  as  a  function  of  temperature  than  simply  assuming  that  k  and  c.  are 
constants.  The  following  will  present  solutions  corresponding  to  the  two  types 
of  boundary  conditions  given  in  Eqs.  (6)  and  (7). 

A.  Prescribed  Heat  Flux  *  at^ 

With  Cjij  «  Cjj  “  1 ,  m  ■  0.04727,  and  n  ■  0.2491  as  given  in  Eqs.  (35)  and 
(36),  the  heat  equation  (25)  subjected  to  the  boundary  conditions  (26)  was 
solved  by  using  the  Adams-Bashf orth  method^  for  b  ■■  0,  1,  2,  and  5.  The 


^R.  Beckett  and  J.  Hurt,  Numerical  Calculations  and  Aleorithms.  McGraw-Hill 
Book  Company,  New  York,  1967,  p.  210. 


Figure  3.  Thermal  Properties  of  Type  420  Stainless  Steel  (normalized) 


second  condition,  f  =  0  at  V  “  ,  was  numerically  satisfied  when  the 

calculation  was  carried  out  to  6.  The  solutions  of  the  function  f  and  its 
derivative  f'  are  tabulated  in  Table  1  and  plotted  in  Figures  A  and  5, 

respectively.  For  comparison,  solutions  of  f  and  f'  for  the  case  of  constant 

properties  (i.e.,  =  Cjj  =  1 ,  m  *  n  =  0)  are  also  provided  in  the  figures. 

Mow  the  local  temperature  0  in  the  solid  can  readily  be  obtained  from  Eq. 
(23).  Figure  6  presents  the  result  for  a  simple  case  of  constant  properties, 
i.e.  K  “  Kq,  C_  •=  CpQ,  and  Qj,  ■  constant.  The  agreement  with  the  exact 
solution  from  Ref.  r  is  excellent.  The  surface  temperature  as  a  function  of 
time  is  represented  by  the  dashed  line  in  Figure  7.  In  another  case  of 
constant  properties,  when  the  boundary  condition  specified  at  the  wall 
surface  is  time  dependent  heat  flux  (Q^  ~  at^),  the  resultant  surface 

temperature  is  shown  in  Figure  8  (dashed  line).  These  figures  show  no  visible 

difference  between  the  present  solutions  obtained  by  the  similarity 
transformation  and  the  exact  solutions.  When  the  temperature  dependence  of 
the  thermal  properties  are  accounted  for  (i.e.,  variable  properties),  the 
exact  solutions  do  not  exist  and  thus  no  comparison  can  be  made.  The  solid 
lines  in  Figures  7  and  8  represent  the  surface  temperature  for  variable 
properties.  They  are  considerably  lower,  especially  at  large  times,  than  the 
surface  temperatures  for  constant  properties.  This  is  due  to  the  fact  that 
heat  is  absorbed  and  conducted  away  faster  when  the  temperature  dependence  of 
the  thermal  properties  is  considered.  The  third  line  in  Figure  7  shows  the 
discrepancy  of  their  results  for  0^  =  constant. 


Table  1.  Solutions  of  Temperature  Function  f  and  Its 
Derivative  f'  at  Surface  (  ’?  »  0) 


Qw  “ 


b 

fc 

fc' 

< 

^v 

0 

1.128A 

-1 .0000 

1.2854 

-0.9882 

1 

0.7523 

-1 .0000 

0.8977 

-1  .0051 

2 

0.6018 

-1.0000 

0.7380 

-1.0145 

5 

0.4168 

-1 .0000 

0.5352 

-1 .0300 

T  c 
w 

T^d  +  t 

b*) 

b* 

fc 

^c' 

fv 

0 

1.000 

-0.5644 

1.000 

-0.6427 

1 

1.000 

-1.1283 

1.000 

-1.0129 

2 

1.000 

-1 .5045 

1.000 

-1.3743 

5 

1.000 

-2.2926 

1.000 

-2.0922 

Subscripts  c  and  v  denote  constant  properties  and 
variable  properties,  respectively. 
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Figure 


Variable  Properties,  Similarity  Solution 
Constant  Properties,  Similarity  Solution 
Constant  Properties,  Exact  Solution 


7.  Surface  Temperature  vs.  Time 

(boundary  condition  Q  ■  constant) 


TIME  (t) 


-  Variable  Properties,  Similarity  Solution 

-  Constant  Properties,  Similarity  Solution 

O  Constant  Properties,  Exact  Solution 


Figure  8.  Surface  Temperature  vs.  Time 
(boundary  condition  Q_  •=  at^) 


B.  Prescribed  Surface  Temperature  T 


Tjjd  +  t°  ) 


With  the  same  values  of  C  ,  C  ,  m,  and  n  used  above,  the  temperature 
function  f  and  its  derivatives  z  solved  from  Eq.  (25)  are  tabulated  in  Table 
1  and  plotted  in  Figures  9  and  10.  Using  those  values,  the  local  temperature 
and  heat  flux  in  the  solid  can  be  calculated  from  £q.  (23)  and  Eq.  (34), 
respectively. 

In  the  case  that  T^  >  constant,  the  temperature  0  should  be  defined  in 
the  form  of  Eq.  (32)  as  explained  earlier.  As  an  example,  if  T^  is  chosen  to 
be  1033°K  (1400*^F),  the  expressions  in  (35)  and  (36)  are  altered  to  the 
following: 

k  -  C  0“  -  1.2060  (37) 

m 

Cp  -  C^0°  -  2.4478  00*328^  (38) 

Since  b*  ~  0  in  this  case,  the  local  temperature  is  simply  given  as 

s 

0  •  f(v)  where  V  » — .  (39) 

Jt 

Figure  11  shows  the  local  temperatures  for  both  constant  and  variable 
properties.  The  temperature  for  constant  properties,  represented  by  the 
dashed  line,  matches  precisely  the  exact  solution.  The  solid  line  next  to  it 
shows  a  lower  temperature  for  variable  properties.  As  explained  in  the 
previous  case,  this  is  due  to  the  fact  that  heat  is  absorbed  and  conducted  to 
the  low  temperature  region  faster  when  the  temperature  dependence  of  the 
thermal  properties  is  taken  into  account.  Figure  11  also  presents  a 
comparison  of  local  heat  fluxes  versus  v  for  the  two  kinds  of  thermal 
properties.  The  quantity  in  the  figure  is  calculated  from 


a  -Q 

“v  Hr  ^Tn*’v 


where  q^  and  q^  are  heat  fluxes  for  constant  properties  and  variable 
properties,  respectively.  Using  the  data  given  in  the  expressions  (37)  and 
(38),  the  calculated  values  of  f^,  f^,  and  are  listed  in  Table  2. 

The  similarity  variable  v,  at  a  given  time,  can  be  treated  as  the  distance 
from  the  surface  of  the  solid.  The  result  in  the  figure  shows  that  at  the 
surface  the  local  heat  flux  q^  for  the  variable  properties  is  greater  than  the 
local  heat  flux  q^  for  the  constant  properties  by  37  percent.  It  is 
interesting  to  note  that  the  percentage  drops  to  zero  at  1  *  1.4  and  then 
turns  negative  as  17  continuously  increases.  To  explain  the  reason  we  note 
that  in  Figure  11  the  two  temperature  curves  merge  at  the  surface  and  also  far 
away  from  the  surface.  Consequently,  as  7  increases  to  a  certain  value,  the 
gradient  of  the  solid  line  changes  from  greater  to  smaller  than  the  gradient 
of  the  dashed  line.  Meanwhile,  the  thermal  conductivity  and  specific  heat  in 
the  expressions  (37)  and  (38)  for  variable  properties  decrease  with 
temperature  when  v  increases.  The  combination  of  these  two  factors  explains 
the  foregoing  change  in  heat  flux. 
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Figure  9.  Profiles  of  Temperature  Function 

(boundary  condition  *  T^d  +  t'’*)) 
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Figure  11.  Temperature  Distribution  and  Heat  Flux 
(boundary  condition  ■  constant) 
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Table  2.  Values  of  Temperature  Function  f  and  Its  Derivatives  f' 
Inside  the  Solid  for  Ty  >  Constant 


V 

fc' 

fv 

^v 

0.00 

1 .00000 

-0.5644 

1.00000 

-0.6427 

0.373 

0.24 

0.86519 

-0.5563 

0.84627 

-0.6333 

0.356 

0.48 

0.73420 

-0.5326 

0.69855 

-0.5926 

0.308 

0.76 

0.59084 

-0.4883 

0.54251 

-0.5172 

0.223 

1.00 

0.47931 

-0.4394 

0.42736 

-0 .4405 

0.138 

1.24 

0.38036 

-0.3841 

0.33094 

-0.3625 

0.052 

1.48 

0.29506 

-0.3263 

0.25269 

-0.2900 

-0.028 

1.76 

0.21302 

-0.2600 

0.18189 

-0.2173 

-0.107 

2.00 

0.15698 

-0.2075 

0.13597 

-0.1664 

-0.161 

2.24 

0.11283 

-0.1609 

0.10100 

-0.1258 

-0.199 

2.48 

0.07915 

-0.1214 

0.07467 

-0.0942 

-0.220 

V.  SUMMARY  AND  CONCLUSIONS 

Similarity  solutions  were  obtained  for  the  transient  heat  conduction  in  a 
semi-infinite  solid  with  temperature-dependent  thermal  properties.  The  method 
of  similarity  transformation  via  one-parameter  groups  provides  a  valuable 
means  for  the  present  analysis.  The  solution  procedure  is  straightforward  and 
the  solutions  obtained  are  highly  accurate,  which  can  be  used  as  reference 
data  for  numerical  solutions  of  similar  nonlinear  problems. 

Results  of  the  temperature  distribution  and  the  local  heat  flux  in  a 
sample  solid  were  obtained  for  two  types  of  time  dependent  boundary 
conditions,  namely  heat  flux  Qy  =  at^  and  temperature  T  =  T^d  +  t^*).  In 
the  limiting  case  of  constant  properties,  the  present  results  are  in  excellent 
agreement  with  existing  exact  solutions.  A  comparison  of  the  results  for 
constant  properties  and  for  variable  properties  has  demonstrated  the 
significance  of  temperature  dependence  of  thermal  properties  for  heat  transfer 
calculations . 
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